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ABSTRACT 


This  paper  examines  the 
problem  of  target  localization 
based  solely  upon  bearing  informa¬ 
tion  obtained  by  a  single  platform 
over  short  observation  times.  A 
number  of  Monte  Carlo  and  analyti¬ 
cal  techniques  for  the  construction 
of  statistical  distributions  of 
target  ranges  are  presented  and 
compared . 


The  actual  range  to  the  target  depends  upon  the  target  motion  and 
the  actual  bearing  errors;  If  the  bearings  are  noisy,  an 
appropriate  question  is:  what  Is  the  distribution  of  range, 
given  the  observed  bearings? 

In  this  paper,  a  number  of  Monte  Carlo  and  analytical 
techniques  for  the  construction  of  the  range  distribution  are 
given.  Sections  2  and  3  contain  descriptions  and  results  of  the 
Monte  Carlo  simulations.  Sections  4  and  5  contain  analytical 
results.  In  section  2  a  “random  range”  simulation  is  described. 
In  section  3,  a  "random  bearing”  simulation  is  described.  In 
section  4,  a  method  for  approximately  calculating  the  mean  and 
variance  of  the  range,  given  the  observations,  is  presented.  In 
section  5,  a  method  for  approximately  calculating  the  distribu¬ 
tion  on  the  range  is  presented.  Section  5  contains  a  discussion 
of  the  tactical  use  of  the  results  presented  in  sections  2  to  4. 

2.  THE  RANDOM  RANGE  SIMULATION 

Description  of  the  Simulation 

This  simulation  runs  for  specified  values  of  the  observed 
initial  and  final  bearings.  The  actual  final  bearings  are 
assumed  to  be  uniformly  distributed  between  +6  and  -6  of  the 


observed  final  bearings.  The  actual  initial  bearings  are  assumed 
to  be  within  ±6  of  the  observed  initial  bearings  with  no 
assumption  as  to  the  distribution.  Final  target  ranges  are 


assumed  to  be  uniformly  distributed  between  Rm^n  and  the  lesser 
of  Rnax  or  a  specified  range  bound  (hence,  the  name  random 
range  simulation.)  Target  motion  is  assumed  to  be  in  a  uniformly 
random  direction  at  a  constant  velocity,  V8« 

The  simulation  begins  with  the  final  observed  bearing.  It 
selects  a  bearing  error  to  determine  the  actual  final  bearing, 
then  selects  a  final  target  range,  R.  It  selects  the  direction 
of  target  motion,  then  traces  motion  of  target  and  observer  back¬ 
wards  in  time  to  their  positions  when  the  initial  bearing  was 
observed.  If  the  actual  bearing  to  the  target  at  this  point  is 
within  +6  or  -6  of  the  observed  initial  bearing,  R  is 
retained  as  a  valid  sample.  Otherwise,  it  is  discarded.  The 
procedure  is  repeated  until  the  desired  sample  size  has  been 
obtained. 

Simulation  Results 

Assume  that  the  observer  is  an  aircraft  and  the  target  is  a 
ship.  The  aircraft  ground  velocity  was  420  kts,  the  target 
velocity  30  kts,  and  the  target  range  was  not  allowed  to  exceed 
75  nautical  miles.  The  time  between  bearing  observations  ranged 
from  30  to  120  seconds  in  15  second  increments.  The  initial 
observed  bearing  was  30  degrees  in  all  instances,  with  bearing 
changes  ranging  from  2  to  10  degrees.  Final  bearing  errors  were 
assumed  to  be  uniformly  distributed  up  to  1  degree  either  side  of 
the  observed  bearing.  One  thousand  replications  were  collected 


for  each  combination  of  bearing  change  and  time  between 
bearings.  Results  of  these  simulations  are  shown  in  and  figures 
2  through  8. 


3.  RANDOM  BEARING  SIMULATION 


Description  of  the  Simulation 


The  simulation  described  in  this  section  works  directly  with 
the  distribution  on  R,  instead  of  assuming  a  distribution  and 
then  rejecting  trials  in  which  the  assumptions  are  not  met. 


Assume  that  the  bearing  changes  are 


and  that  the 


aircraft  and  ship  speeds  are  VA  and  V§,  that  the  observation 
time  is  t  ,  and  that  the  target  moves  in  the  direction  a.  In 
the  next  section  it  is  shown  that  the  range  R  to  the  target 
after  the  second  bearing  is 


VAt  sin  (4»q) 
sin  -  ♦q) 


Vgt  sin  (<|>0  -  a) 
sin  -  <J»0) 


(3.1) 


The  simulation  proceeds  as  follows: 

A  A 

1.  Initialize,  N-0.  Input  observed  bearings  <b  q  » 


2.  Select  a,  where  a  is  uniformly  distributed  on  [0,  360°] 


TABLE  3 


FRACTION  OF  TRIALS  WITH  RANGE  GREATER  THAN  R* 


420  kts  Vg  -  30  kta  t  ■  60  seconds  R  <  75  n.nl 


Bearing  change  (degrees 


TABLE  S 


R  =  60 


RMAX  <  75  n.mi. 
rMAX<  150n.mi. 


Observed  bearing  change  (degrees) 


FIG.  2:  DISTRIBUTION  OF  TARGET  RANGES  FOR  OBSERVED  BEARING  CHANGES 
WITH  VA  -  420  KTS,  V,  -  30  KTS,  t  -  30  SECONDS,  AND  BEARING  ERRORS 
UNIFORMLY  DISTRIBUTED  BETWEEN  ±1°  AND  RMAX  <  75  n.mi.  (1,000 
REPLICATIONS  PER  BEARING  CHANGE  WITH  INITIAL  BEARING  OF  30°) 


"MAX  <  75  n-ml 


Observed  bearing  change  (degrees) 


FIG.  3:  DISTRIBUTION  OF  TARGET  RANGES  FOR  OBSERVED  BEARING  CHANGES 
WITH  VA  -  420  KTS,  V,  -  30  KTS,  t  -  45  SECONDS,  AND  BEARING  ERRORS 
UNIFORMLY  DISTRIBUTED  BETWEEN  ±1°  AND  RMAX  <  75  n.mi.  (1,000 
REPLICATIONS  PER  BEARING  CHANGE  WITH  INITIAL  BEARING  OF  30°) 


Observed  bearing  change  (degrees) 


FIG.  5:  DISTRIBUTION  OF  TARGET  RANGES  FOR  OBSERVED  BEARING  CHANGES 
WITH  VA  -  420  KTS,  V,  -  30  KTS,  t  -  75  SECONDS,  AND  BEARING  ERRORS 
UNIFORMLY  DISTRIBUTED  BETWEEN  ±1°  AND  RMAX  <  75  n.mi.  (1,000 
REPLICATIONS  PER  BEARING  CHANGE  WITH  INITIAL  BEARING  OF  30°) 


Observed  bearing  change  (degrees) 


FIG.  7:  DISTRIBUTION  OF  TARGET  RANGES  FOR  OBSERVED  BEARING  CHANGES 
WITH  VA  -  420  KTS.  V,  -  30  KTS,  t  -  105  SECONDS,  AND  BEARING  ERRORS 
UNIFORMLY  DISTRIBUTED  BETWEEN  ±1°  AND  RMAX  <  75  n.mi.  (1,000 
REPLICATIONS  PER  BEARING  CHANGE  WITH  INITIAL  BEARING  OF  30°) 


Observed  bearing  change  (degrees) 


FIG.  8:  DISTRIBUTION  OF  TARGET  RANGES  FOR  OBSERVED  BEARING  CHANGES 
WITH  VA  -  420  KTS,  V,  -  30  KTS,  t  -  120  SECONDS,  AND  BEARING  ERRORS 
UNIFORMLY  DISTRIBUTED  BETWEEN  11°  AND  RMAX  <  75  n.mi.  (1,000 
REPLICATIONS  PER  BEARING  CHANGE  WITH  INITIAL  BEARING  OF  30°) 


is  computed  directly.  In  the  random  range  simulation,  the 
distribution  function  on  R  is  Initially  assmed  to  be  uniform 
on  timin'  *«axJ  (where  Rmax  is  either  Rmax  from  equation 
(1.2)  or  an  artificially  Imposed  cut-off).  Then  trials  that  do 

A 

not  conform  to  the  assumption  that  -1<  $0  ~  $o  are 

rejected.  The  differences  in  the  two  procedures  should  be 
reflected  most  in  differences  in  the  tails  of  the  range 
distributions. 


Figure  9  shows  the  results  of  the  two  kinds  of  simulations 
It  appears  that  for  most  ranges,  the  R  simulation  underesti¬ 
mates  the  closeness  of  the  target. 


4.  ANALYTICAL  RESULTS:  MEAN  AND  VARIANCE  OF  THE  RANGE 


First  consider  the  derivation  of  equation  (3.1).  The 


geometry  is  shown  in  figure  10.  Set  S  ■  VAt  and  d  ■  Vgt 
Applying  the  law  of  sines  to  the  small  triangle  gives 


For  the  large  triangle,  the  law  of  sines  gives 


s in ( ) 


Solving  equations  (4.2)  and  (4.3)  for  R  gives 


I 

S  sin(<frQ) 

R  - - Q 

sin(<|)1  -  4>0) 

I 

S  sin($n)  d  sin($n-a) 

-  - - - - - -  (4 

sin(«t>1  -  4»0)  sin(^  -  <J»0) 

I 

Now  replace  $q  and  ^  with  random  variables,  X  and  Y. 
is  a  random  variable  too.  Furthermore,  R  is  a  nonlinear 
I  function  of  the  random  variables  X  and  Y.  In  order  to 

calculate  the  mean  and  variance  of  R,  the  following  approxi¬ 
mation  is  used.  If  f (X ,Y)  is  any  twice  differentiable  func 
[  tion,  and  X  -  E{X>,  X  -  E{Y} ,  then  a  Taylor  expansion  gives 

the  approximation 


E{f(X,Y)}  IB  E{f(X,Y)  +  f  (X,y)  { X-X]  +  f  (X  ,Y)  [Y-Y] 

X  / 


+  Y  f xx (X , Y )  [X-Y]2  +  fxy(X,Y)  [X-X11Y-Y]  (4 

I 

+  i  f  „ „ (  X  ,  Y )  [  Y-Y  ]  2 } 


so  that 


E{  f  (X  ,  Y)  }  HI  f(X,Y)  +  j[fxx(X,Y)  Var  X 


(4.5) 


+  2f  (X , Y)  Cov(X.Y)  +  f  (X,Y)  Var  Y- 

xy  y  y 

In  the  same  way  one  finds  that 

E{f2(X,Y)}  I  f 2(X,Y)  +  { f 2 (X , Y)  +  f (X  ,  Y)  fxx(X,Y)))  Var  X 

+  [fj  (X,Y)  +  f(X,Y)  f yy (X ,  Y  )  ]  Var  Y 

(4.6) 

+  2  (  f  x ( X  ,  Y)  fy(X,Y) 

+  f  (  X  ,  Y  )  f  (X,Y)]  Cov  (X  ,  Y ) 
xy 

In  the  case  of  interest  here,  f(X,Y)*R“[S  sin(X) 

-  d  sin(X-a) ] [ s in( Y-X) ] “ *  ,  and  assuming  Cov(X,Y)*0,  the 
derivatives  of  interest  are 


-25- 


Ry  -  -[S  sln(x)  -  d  sin(x-a)*  [sin(y-x)»  cos(y-x) 

■  -  [s  sin(x)-d  sin(x-a)"  {-2[sin  (y-x)»  cos4(y-x) 
-  [sin  (y-x)-  1 } 


( 


Rx  »  [S  cos (x)-d  cos(x-a)«  [sin(y-x)* 

_  2 

+  [s  sin(x)  -  d  sin(x-a)"  [sin(y-x)-  cos  (y-x) 

.  2 

Rxx  ■  2tS  cos(x)  -  d  cos(x-a)*  [sin(y-x)»  cos(y-x) 

”3  2 

+  2 [s  sin(x)  -  d  sin(x-a)*  [sin(y-x)-  cos  (y-x) 


The  derivatives  (4.7)  involve  cos(x-o)  and  sin(x-o). 
Assuming  the  a~U[wg,wjJ,  table  1  shows  the  averages  of  the 


functions  of  interest. 


cos(x-a) 

0 

360° 

0 

s i n ( x-a) 

0 

360* 

0 

cos (x-a) 

0 

360* 

2  sinx/x 

s in(x-a) 

0 

360* 

-2  cosx/x 

Note  that  If  a  is  uniformly  distributed  on  [0,  360*1,  then  the 
target  motion  makes  essentially  no  contribution  to  the  results  of 
this  level  of  approximation.  Using  equations  (4.5)  -  (4.7)  and 
table  8,  assuming  X  and  Y  are  independent,  and  X,Y  ~  Ul~l*, 

1°1  leads  to  figu  res  11,12,  and  13.  In  these  figures,  the  mean 
of  R  is  plotted  against  f°r  *  “  30*  40,  and  50 

seconds,  with  a  ~  U[G,2n]»  ■  420  kts ,  and  Vg  ■  30  kts.  The 

error  bars  indicate  one  standard  deviation. 

The  considerable  size  of  the  error  bars  suggests  that 
distributional  results  are  very  important,  since  knowledge  of  the 
mean  and  variance  is  not  sufficient  in  this  case  to  give  accurate 


predictions 


5.  ANALYTICAL  RESULTS:  APPROXIMATE  CALCULATION  OF  THE 
DISTRIBUTION  ON  R 

Since  R  Is  a  nonlinear  function  of  X  and  Y,  It  Is 
difficult  to  calculate  its  distribution  exactly.  In  this 
section,  an  approximate  method  for  the  calculation  of  the 
distribution  is  given.  The  starting  point  is  equation  (4.3) 
for  R  with  X  and  Y  replacing  an^  4i: 

S  sln(X)  d  sin(X-a) 

R  - -  (5.1 

sin(Y-X)  sin(Y-X) 

The  approximation  used  here  is  the  small  noise  approximation,  ii 
which  X  and  Y  are  replaced  by 
X  -  X  +  6 

(5.2] 

y  -  y  +  n 

where  5~U[ -6^  ,  61  ]  and  q~U  I ,t]^  ] 

are  the  deviations  from  the  average,  and  a  Taylor  expansion  of 
(5.1)  is  used  to  write  R  as  a  linear  function  of  6  and  q. 
Such  a  procedure  is  valid  if  5  and  T)  are  sufficiently  small 
relative  to  other  terms  in  the  expansion.  Although  used  mainly 
for  mathematical  convenience,  if  such  an  approximation  is  not 
valid  the  entire  concept  of  getting  useful  information  from  two 
noisy  bearings  may  not  be  valid,  since  then  the  noise  overrides 


the  deterministic  terms 


Error  bars  indicate  1 
standard  deviation 


In  this  approximation 


sin(X)  1®  8  in  ( X )  +  6  cos(X)  =  +  Aj6 

sin(X-a)  Ht  sin(X-a)  +  6  cos(X-a)  =  +  C25 

[sin(Y-X)]'1  IB  [sin(Y-X)  ]_1 

cos ( Y-X) 

- 1 — I - r  (tj-6)  s  B.  -B„  (n-6) 

[sin(Y-X)]2  1  2 

Using  (5.3)  in  (5.1)  gives 

R  -  R  +  E6  +  Ft) 

where 

R  -  B1(SA1  -  dCl) 

E  -  B,(SA,  -  dC_  )  +  B«  ( SA,  -  dC,  ) 


The  distributional  quantity  of  interest  H(r)  is  then 


H(r )  «  Pr{R  <  r}  ■  PrfR  +  E6  +  Ft)  <  r} 


(5.6 


It  is  easy  to  show  that,  for  the  operational  values 
interest,  E>0  and  F<0.  In  light  of  the  form  of 
there  are  maximum  and  minimum  possible  values  of  R 


of 

R  in  (5.4) 
given  by 


R  -  R  +  E6.  +  IfU, 
max  1  1  1  1 


min 


R  -  E61  +  Fri] 


(5. 


Note  that  Rmax  and  ^min  depend  upon  a  through  the 

coefficients  In  (5.3)  and  (5.5). 


Define  new  random  variables,  U  and  V  by 

U  -  E6  V  ■  |  F  1 1]  (5.8 


so  that 


In  the  sequel  it  is  helpful  to  set  u^  ■  E6j,,  v^  »  i  ^  I  'H 1 
Defining  Z  5  r-R,  equation  (5.6)  becomes 


Pr { R  <  r}  -  Prfu  -  V  <  z} 


(5.10) 


In  order  to  compute  (5.10),  refer  to  figure  14.  For  a 
given  Z,  Fr{U  -  V  <  z}  is  the  fraction  of  the  area  of  the 
rectangle  above  the  line  u  -  v  »  Z.  There  are  three  cases.  Set 
Zi  -  ui  -  Vl  and  z2  “  V1  “  ul»  then  the  three  cases  are 
Case  I  :  Z> Zj 
Case  II:  Z<Z2 
Case  III:  Z2<Z<ZL 


In  Case  I 


Pr {u  -  V  <  z}  -  1  - 


8ulVl 


(U1  +  Vj  -  Z) 


(5.11) 


In  Case  II, 


Pr {u  -  V  <  Z}  - 


8ul  vl 


(Z  +  ux  +  vL)' 


(5.12) 


In  Case  III, 


Pr {u  -  V  <  Z}  -  1  - 


{ 2v ?  +  2v . ( u. -  Z  -  v. )}  (5.13) 


Figure  15  shows  a  comparison  of  the  theory  developed  here 
with  the  range  and  bearing  simulations  for  the  same  parameters  as 
in  figure  9.  Further  comparisons  of  the  theory  and  bearing 
simulation,  with  the  noise  uniformly  distributed  on  dfc  .5°  or  ± 
•25°  showed  that  the  agreement  between  theory  and  simulations 
was  very  good. 

Further  comparisons  of  the  theory  and  the  random  range 
simulation  showed  that  the  predicted  minimum  values  for  the 
ranges  agreed  closely,  but  that  the  predicted  maximum  values 
disagreed.  As  an  example,  for  the  parameter  values  in  figure  9, 
the  minimum  range  by  the  theory  is  34  n.ml.  and  the  minimum  range 
by  the  range  simulation  is  35  n.mi.  The  maximum  ranges,  however, 
are  86  n.mi.  and  117  n.mi.  respectively.  The  cause  of  the  long 
tail  on  the  range  simulation  is  not  known. 

TACTICAL  IMPLICATIONS  AND  USES 

One  tactical  use  of  the  techniques  described  here  is  the 
determination  of  confidence  levels  on  the  range  to  a  target. 
Referring  to  figure  16,  if  the  observer  wished  to  know  that  he 
was  within  60  miles  of  the  target  with  .9  probability,  he  would 
need  to  observe  a  bearing  after  t  seconds  that  was  above  the 
Pgg  -  .9  curve.  If  the  observer  wished  to  know  that  he  was 
outside  30  miles  of  the  target  with  .9  probability,  he  would 
need  to  observe  a  bearing  after  t  seconds  that  was  below  the 
P30  ■  .1  curve.  Finally,  if  the  observer  wished  to  know  that  he 


t,  time  from  initial  bearing  (seconds) 


FIG.  16:  PA0  AND  P30  PROFILES  FOR  BEARING  OBSERVED  t  SECONDS 
AFTER  THE  INITIAL  BEARING  OF  20  DEGREES 


was  somewhere  between  30  and  60  miles  from  the  target  with  .9 
probability,  his  bearing  observation  after  t  seconds  would  need 
to  be  below  the  Pjq  ■  .1  curve  but  above  the  Pgg  “  *9 
curve.  This  is  the  shaded  region  shown  in  the  figure. 

There  may  be  other  applications  for  this  sort  of  informa¬ 
tion.  For  any  application  to  be  operationally  useful  it  will  be 
necessary  to  determine  through  operational  testing  the  actual 
distribution  of  bearing  errors  encountered  in  the  bearing 


observations 
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